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Abstract 

The origin of fluctuations in the average number of intermediate mass fragments seen in experi- 
ments in smah projectile hke fragments is discussed. We argue that these can be explained on the 
basis of a recently proposed model of projectile fragmentation. 
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1 



This report is an outgrowth of a recent work |l| where we introduced a model of projectile 
fragmentation in heavy ion collisions. The model has three ansatzs: at a given impact 
parameter a certain part of the projectile is sheared off forming a projectile like fragment 
(PLF) with charge number and neutron number A^^. This is called abrasion. This hot PLF 
expands to about one-third the normal density and then breaks up into several composites 
at a given temperature T. This break up is according to a canonical thermodynamical model 
(CTM). The composites are hot and will decay by evaporation. For details, please see 
The main contention of reference [l| was that the temperature T must be taken to be a 
function of the impact parameter h. 

Results of the following experiment done at the SIS heavy-ion synchrotron at GSI 
Darmstadt are published |2|. In an event let us denote the number of intermediate mass 
fragments (IMF) (charge z between 3 and 20) by Nimf- In the same event denote by 
■2^f)onnd=sum of all the charges in the PLF minus the charges oi z = 1 particles (proton, 
deuteron and triton). After many events one can plot Mjmf vs Zbound where Mjmf is the 
average of Njmf- The data and comparison with the theoretical calculation done in jl| is 
shown in Fig. 1. 

The overall feature of the figure is that the general shapes of the theoretical and 
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TABLE L Experimentally measured Mjmf at small Z^ound values for Sn^^''^ on Sn^^^, Sn^'^^ on 
Sn^^^ and La^'^^ on Sn^^'^ reactions. 



experimental curves agree. The b dependence of T is crucial for this (as explained in 
[l|). However, there are significant fluctuations in the experimental values of Mjmf for 
low values of Zf,ound whereas theory completely misses these fluctuations. In this note we 
explain (a)why these fluctuations arise, (b) how staying within the main ingredients of the 



theoretical model but using more realistic parameters we can reproduce the fluctuations 
and (c) why the calculation in l| missed the fluctuations seen in small PLF's. 



First we explain how the fluctuations arise. We have Zhound=Zs minus the sum of charges 
of all z=l particles (protons, deuterons and tritons). Also a particle is considered to be 
an IMF if its charge z is between 3 and 20. Just these two conditions and some general 
knowledge of low-mass nuclei allow us to reach some interesting conclusions. 

If Zbound='i it guarantees that we have a 3 Li nucleus. Thus for Zbound='i Mjmf=^- If 
3 Li decays by a proton emission we are no longer in Zhound='i but degenerate into Zhound='^- 
Also there is no IMF. If it decays by neutron emission to a particle stable state of a different 
isotope of Li, we still have Mimf=^- There are several particle stable states of Li so 
Zbound = 3, M/MF=1 is always satisfled. 

Let us consider now Z}yound=^- For Zbound=^ we can have a Be nucleus with Njmf=^ but 
it can also decay into two He isotopes which still retains Zbound=4: but with Nimf=0- We 
therefore expect to have Zbound=4: and Mjmf = X where X is less than 1. But unfortunately 
X is not the same value for all Be nuclei. We will soon demonstrate how X could be 
determined for each Be nucleus but the fact that X varies from one isotope of Be to another 
isotope of Be makes the evaluation of Mjmf in the case of Zbound=4: a lengthy procedure. 

If Zbound=5 we have either one Boron nucleus or a Li nucleus plus a He nucleus. In both 
the cases Mimf=^- If the Boron nucleus sheds a proton, the status drops to Zf)ound=4: and 
we are back to the Zbound=4: case. If the Boron nucleus sheds one or more neutrons to reach 
a particle stable state we maintain Zbound = 5, Mimf=^- If Boron decays into a Li and 
He two things can happen. We reach a particle stable state of Li and we have Zbound = 5, 
MiMF=^- If the Li sheds a proton we no longer have .^60^^=5. Thus so long as we have 
Zbound=5 we have M/mf=1- 

We want to get back to the case of Zbound=^- Now we need to bring in details of the 
model. Two modiflcations are made. To carry out CTM one needs to put in the partition 
function of each composite into which the hot abraded PLF can break into. In our previous 
calculation, except for nuclei upto '^He, we used the liquid-drop model for the ground state 
energy and the Fermi-gas model for excited states. For small PLF's this is inaccurate and 
we put in experimental values of ground state and excited state energies. Usually all excited 
states upto 7.5 MeV are included. Next we consider the decays of hot composites resulting 
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from CTM. Previously we used an evaporation code. We replace this by actual decay data 
whenever possible. In practical terms this means the following. A nucleus has many energy 
levels and a hot nucleus means that the probability of occupation of a state i is proportional 
to Siexp{—exc{i) /T)) where Si is the spin degeneracy and exc{i) is the excitation energy. 
The decay of the state i is taken from data table l3-7| where available or guessed from 

n 

systematics. We take T to be 7 MeV suggested by our past work 

It is useful to list first the deacy properties of hot Be nuclei. These are computed at T=7 
MeV. 

n 

^Be: This decays into [3] gHe plus 2 protons so this counts as Zhound='^ and Mjmf=0- 



^Be: The lowest 2 states are particle stable. Population into any of these gives Zi 



bound 



and one IMF. The probability of this occurring is 0.406. The other states decay into pie plus 
gHe leading to Zbound=4: and Nimf=0. Thus for jBe we have Zf,ound=4: and M/mf=0.406. 

4Be: this occurs as resonances of two 2He so here Z]jound=^ and Mimf=^- 

4Be: Only the ground state is particle stable, the rest decay to neutron plus two alphas. 
The occupation probability in the ground state is 0.193. So here Z]jound=^ and M/AfF=0.193. 

4*^Be: Here we have taken all the levels upto 6.26 MeV (summed occupation probabil- 
ity=0.604) to give Zi,ound=^ and 1 IMF and rest of the levels upto 9.3 MeV to give Zbound=^ 
and IMF. Thus ^^^=4 and M/mf=0.604. 

4^Be: Here the lowest two levels have Zbound=4: and Njmf=^ and the probability of 
occupation 0.1567. The higher levels, with summed occupation probability 0.8433 go to 
4°Be+n. We have assigned them Zbound=^ and M/jv/i?=0.604. Thus we take 4^Be to give 
2'6o«nd=4 and M7-mf=0.666. 

Let us now outline how we calculate Mjmf for Zbound=4: for collisions of and 
^^^La on ^^^Sn. Although our discussion will be limited to Zf,ou„d=4, the method can be 
extended to higher values of Zhound except that the complexity increases very rapidly. The 
method of obtaining the abrasion cross-section for a PLF with given Zs,Ns is given in [l|. 
For Zhound=4: we need to consider Zs=A (most important) and higher. Once a PLF with 
given Zg, Ng is formed it will expand to one-third the normal nuclear density and break up 
into hot composites. Just as we could characterize a hot Be nucleus by a Zbound and Mjmf 
we can ascribe to each Zs,Ns a probability of obtaining Zbound=^ with an associated Mjmf- 
(An example below shows how this can be done.) Table II compiles these values (last two 
columns) . 
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TABLE II: Abrasion cross-sections (in millibarns) for a given (Zg, Ng) for Sn , Sn^'^^ and La 
on Sn^^^. P{Zi,ound = 4) gives the probability of obtaining Zio^nd = 4 for a given Zg, Ng and 
MiMF is the corresponding average multiplicity of intermediate mass fragments. 

Utilizing also the values of the abrasion cross-sections for Zg, Ng for the three reactions (also 
given in the Table) we get the desired results. For Z^ound—^, M/mf=0. 145(0. 14) for -'^^'''Sn 
beam, 0.151(0.178) for -"^^^La beam and 0.38(0) for -"^^^Sn beam. The experimental values are 
enclosed by parenthesis. Except for ^^^Sn beam our results approximately correspond to the 
experimental data. The value for ^^^Sn is a mystery. In any model we can think of the 
result should not be or very different from the other two. In any case we have reproduced 
the fluctuation: Mimf drops from 1 at Zhound—^ to much lower value at Z^ound—^ and back 
again to 1 at Zbound—^- It is very long to do a quantitative estimate for Zhound—^- This 
will arise from Zs=6 and higher. A study of the CTM results of Zg = 6, Ng = 7 shows the 
following. There is a significant probability of reaching a Carbon nucleus {Zf,ound=Q) ■ This 
will produce a Mimf ~1- There is a comparable probability of obtaining Zbound=Q with a 
^Be nucleus (zero IMF) and another He nucleus and also a ^Be nucleus (M7mf=0.193) and 
another He nucleus. The probabihty of reaching two Li nuclei post CTM is non-neghgible 
but the chances of any one or both of them decaying by alpha or proton emission (thereby 
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dropping below Zbound—(i) are quite high (0.88). A theoretical value for Mjmf ~0.5 seems 
quite plausible. It is the fragility of Be nuclei which produces the dip in Mimf for Z^jound—^ 
and is also responsible for the dip at Z(,o„„d=6. 
As promised, let us give an example how for a given Z^, A^^ the probability of occurrence 
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TABLE III: Multiplicity of different fragments produced by CTM from the abraded nucleus = 4, 
A, = 5 at T = 7.0 MeV. 

of Zbound—^ and the associated Mjmf can be computed (last two columns of Table II). 
Consider = 4, A^^ = 5. To start with, the average numbers of each composite resulting 
from the CTM break up of = 4, A^^ = 5 system are listed in Table III. But in a simple case 
like this, this can also give, with little effort, the probability of a channel or the probability 
of a sum of channels. From the table, the average number of 4Be is ^0.475. This is a 
channel where only 4Be and nothing else appears. Thus there is a probability of 0.475 of 
reaching Z^cmnd—^ and M/mf=0.193. Next, looking at the table, the average number of 
^Be is ~0.088. This comes from a channel where there is one ^Be and one neutron. So we 
have a probabihty of 0.088 of reaching Zj,ound—4^ with M/mf=0.0. Next from the table, the 
average number of fLi is ~0.038. This has to occur in combination with a proton. Clearly, 
this is channel with Zhound=^, so this does not concern us presently. Next, from the table, 
the average number of jBe is w0.006. This is a channel which has one jBe and 2 neutrons. 
Thus we have a probability of 0.006 of reaching Z^ound—^ with Mjmf— 0.406. We have 
exhausted all the channels for reaching Zi,ound—^- Summing up with appropriate weightage, 
from Zg — Ng — 5 the probability of reaching .^bound— 4 is 0.569 with M/MF=0.165. 
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We can repeat similar arguments for other Zs, Ng in Table II. The cases of Zs=5 are more 
complicated. 

We now try to answer why the calculation of [l| failed to produce any fluctuation. There 
are many reasons (use of liquid-drop model and non-recognition of the fragility of Be nucleus 
etc.) but the most interesting reason is different. 

The prescription we used for Mjmf vs. Zbound is the following. At a given b, abrasion gives 
an integral Zg (and an integral A^^). This system expands, then dissociates by CTM and the 
hot composites which are the end results of the CTM, can evaporate light particles to give 
the final distribution. From this we obtained Mjmf and we considered Zhound to be given by 
Zg — J2i nz=i{i) where 71^=1 (i) stands for the average multiplicity of proton/deuteron/triton. 
This prescription does not match exactly the experimental procedure. Experimentally Zj^ound 
is obtained event by event and in every event Zbound is an integer (sum of all charges from 
PLF minus number of particles with z=l). From many events with the same Zbound one can 
obtain Mjmf- In our calculations although Zg is an integer Zg — J2i^z=i{i) will usually be 
non-integer since the n2=i(z)'s (average number of composite i) are. 

Our calculation can map much better into a different experiment. In this experiment 
Zg is measured but z=l particles are not subtracted. One then obtains Mjmf for each 
Zg. This problem is simpler: given a total number of particles, what is Mjmf"^ But in the 
reported experiment one asks a more exclusive question : when the particles are fractured 
in a certain way (a given number of particles with charge greater than 1) what is M/^v/f? 
In our prescription we get a non-integral value for Zbound and what we are obtaining is an 
average of Mjmf done over Mjmf belonging to different but neighbouring values of integral 
Zbound- This would be quite wrong if values of Mjmf belonging to neighbouring Zbound s 
differ strongly (as it happens for very small systems) but for large systems the difference 
would be small and our prescription is adequate for an estimate. 
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FIG. 1: (Color Online) Mean multiplicity of intermediate-mass fragments Mjmf, as a function 
of Ziiound for (a) ^^''Sn on ^""^^Sn and (b) ""^^^Sn on ^""^^Sn reaction (red solid lines). Temperature is 
impact parameter (5) dependent and falls off linearly with b from 7.5 MeV at 6=0 to 3 MeV at 



maximum value of b. The experimental results 



are shown by the black dashed lines. 
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